The influence of tree genus, phylogeny, and richness on the specificity, rarity, and diversity of ectomycorrhizal fungi.

Partner specificity is a well-documented phenomenon in biotic interactions, yet the factors that determine specificity in plant-fungal associations remain largely unknown. By utilizing composite soil samples, we identified the predictors that drive partner specificity in both plants and fungi, with a particular focus on ectomycorrhizal associations. Fungal guilds exhibited significant differences in overall partner preference and avoidance, richness, and specificity to specific tree genera. The highest level of specificity was observed in root endophytic and ectomycorrhizal associations, while the lowest was found in arbuscular mycorrhizal associations. The majority of ectomycorrhizal fungal species showed a preference for one of their partner trees, primarily at the plant genus level. Specialist ectomycorrhizal fungi were dominant in belowground communities in terms of species richness and relative abundance. Moreover, all tree genera (and occasionally species) demonstrated a preference for certain fungal groups. Partner specificity was not related to the rarity of fungi or plants or environmental conditions, except for soil pH. Depending on the partner tree genus, specific fungi became more prevalent and relatively more abundant with increasing stand age, tree dominance, and soil pH conditions optimal for the partner tree genus. The richness of partner tree species and increased evenness of ectomycorrhizal fungi in multi-host communities enhanced the species richness of ectomycorrhizal fungi. However, it was primarily the partner-generalist fungi that contributed to the high diversity of ectomycorrhizal fungi in mixed forests.

DNA-metabarcoding supports trophic flexibility and reveals new prey species for the Galapagos sea lion.

Tropical ecosystems are challenging for pinnipeds due to fluctuating food availability. According to previous research, the Galapagos sea lion (GSL, Zalophus wollebaeki) adopts trophic flexibility to face such conditions. However, this hypothesis comes from studies using traditional methods (hard-parts analysis of scat and isotopic analysis from tissue). We studied the diet of five rookeries in the southeastern Galapagos bioregion (which harbors the highest GSL density), via DNA-metabarcoding of scat samples. The DNA-metabarcoding approach may identify consumed prey with a higher taxonomic resolution than isotopic analysis, while not depending on hard-parts remaining through digestion. Our study included five different rookeries to look for evidence of trophic flexibility at the bioregional level. We detected 98 prey OTUs (124 scats), mostly assigned to bony-fish taxa; we identified novel prey items, including a shark, rays, and several deep-sea fish. Our data supported the trophic flexibility of GSL throughout the studied bioregion since different individuals from the same rookery consumed prey coming from different habitats and trophic levels. Significant diet differentiations were found among rookeries, particularly between Punta Pitt and Santa Fe. Punta Pitt rookery, with a more pronounced bathymetry and lower productivity, was distinguished by a high trophic level and consumption of a high proportion of deep-sea prey; meanwhile, Santa Fe, located in more productive, shallow waters over the shelf, consumed a high proportion of epipelagic planktivorous fish. Geographic location and heterogeneous bathymetry of El Malecon, Española, and Floreana rookeries would allow the animals therein to access both, epipelagic prey over the shelf, and deep-sea prey out of the shelf; this would lead to a higher prey richness and diet variability there. These findings provide evidence of GSL adopting a trophic flexibility to tune their diets to different ecological contexts. This strategy would be crucial for this endangered species to overcome the challenges faced in a habitat with fluctuating foraging conditions.

Dataset from RNAseq analysis of differential gene expression among developmental stages of two non-marine ostracodes.

We contribute transcriptomic data for two species of Ostracoda, an early-diverged group of small-sized pancrustaceans. Data include new reference transcriptomes for two asexual non-marine species (Dolerocypris sinensis and Heterocypris aff. salina), as well as single-specimen transcriptomic data that served to analyse gene expression across four developmental stages in D. sinensis. Data are evaluated by computing gene expression profiles of the different developmental stages which consistently placed eggs and small larvae (at the stage of instar A-8) similar to each other, and apart from adults which were distinct from all other developmental stages but closest to large larvae (instar A-4). We further evaluated the transcriptomic data with two newly sequenced low-coverage genomes of the target species. The new data thus document the feasibility of obtaining reliable transcriptomic data from single specimens - even eggs - of these small metazoans.

Connecting the multiple dimensions of global soil fungal diversity.

How the multiple facets of soil fungal diversity vary worldwide remains virtually unknown, hindering the management of this essential species-rich group. By sequencing high-resolution DNA markers in over 4000 topsoil samples from natural and human-altered ecosystems across all continents, we illustrate the distributions and drivers of different levels of taxonomic and phylogenetic diversity of fungi and their ecological groups. We show the impact of precipitation and temperature interactions on local fungal species richness (alpha diversity) across different climates. Our findings reveal how temperature drives fungal compositional turnover (beta diversity) and phylogenetic diversity, linking them with regional species richness (gamma diversity). We integrate fungi into the principles of global biodiversity distribution and present detailed maps for biodiversity conservation and modeling of global ecological processes.

Structure of plant-associated microeukaryotes in roots and leaves of aquatic and terrestrial plants revealed by blocking peptide-nucleic acid (PNA) amplification.

Studies of plant-microbe interactions, including mutualistic, antagonistic, parasitic, or commensal microbes, have greatly benefited our understanding of ecosystem functioning. New molecular identification tools have increasingly revealed the association patterns between microorganisms and plants. Here, we integrated long-read PacBio single-molecule sequencing technology with a blocking protein-nucleic acid (PNA) approach to minimise plant amplicons in a survey of plant-eukaryotic microbe relationships in roots and leaves of different aquatic and terrestrial plants to determine patterns of organ, host, and habitat preferences. The PNA approach reduced the samples' relative amounts of plant reads and did not distort the fungal and other microeukaryotic composition. Our analyses revealed that the eukaryotic microbiomes associated with leaves and roots of aquatic plants exhibit a much larger proportion of non-fungal microorganisms than terrestrial plants, and leaf and root microbiomes are similar. Terrestrial plants had much stronger differentiation of leaf and root microbiomes and stronger partner specificity than aquatic plants.

A pile of pipelines: An overview of the bioinformatics software for metabarcoding data analyses.

Environmental DNA (eDNA) metabarcoding has gained growing attention as a strategy for monitoring biodiversity in ecology. However, taxa identifications produced through metabarcoding require sophisticated processing of high-throughput sequencing data from taxonomically informative DNA barcodes. Various sets of universal and taxon-specific primers have been developed, extending the usability of metabarcoding across archaea, bacteria and eukaryotes. Accordingly, a multitude of metabarcoding data analysis tools and pipelines have also been developed. Often, several developed workflows are designed to process the same amplicon sequencing data, making it somewhat puzzling to choose one among the plethora of existing pipelines. However, each pipeline has its own specific philosophy, strengths and limitations, which should be considered depending on the aims of any specific study, as well as the bioinformatics expertise of the user. In this review, we outline the input data requirements, supported operating systems and particular attributes of thirty-two amplicon processing pipelines with the goal of helping users to select a pipeline for their metabarcoding projects.

Metabarcoding of soil environmental DNA to estimate plant diversity globally.

Introduction: Traditional approaches to collecting large-scale biodiversity data pose huge logistical and technical challenges. We aimed to assess how a comparatively simple method based on sequencing environmental DNA (eDNA) characterises global variation in plant diversity and community composition compared with data derived from traditional plant inventory methods. Methods: We sequenced a short fragment (P6 loop) of the chloroplast trnL intron from from 325 globally distributed soil samples and compared estimates of diversity and composition with those derived from traditional sources based on empirical (GBIF) or extrapolated plant distribution and diversity data. Results: Large-scale plant diversity and community composition patterns revealed by sequencing eDNA were broadly in accordance with those derived from traditional sources. The success of the eDNA taxonomy assignment, and the overlap of taxon lists between eDNA and GBIF, was greatest at moderate to high latitudes of the northern hemisphere. On average, around half (mean: 51.5% SD 17.6) of local GBIF records were represented in eDNA databases at the species level, depending on the geographic region. Discussion: eDNA trnL gene sequencing data accurately represent global patterns in plant diversity and composition and thus can provide a basis for large-scale vegetation studies. Important experimental considerations for plant eDNA studies include using a sampling volume and design to maximise the number of taxa detected and optimising the sequencing depth. However, increasing the coverage of reference sequence databases would yield the most significant improvements in the accuracy of taxonomic assignments made using the P6 loop of the trnL region.

How, not if, is the question mycologists should be asking about DNA-based typification.

Fungal metabarcoding of substrates such as soil, wood, and water is uncovering an unprecedented number of fungal species that do not seem to produce tangible morphological structures and that defy our best attempts at cultivation, thus falling outside the scope of the International Code of Nomenclature for algae, fungi, and plants. The present study uses the new, ninth release of the species hypotheses of the UNITE database to show that species discovery through environmental sequencing vastly outpaces traditional, Sanger sequencing-based efforts in a strongly increasing trend over the last five years. Our findings challenge the present stance of some in the mycological community - that the current situation is satisfactory and that no change is needed to "the code" - and suggest that we should be discussing not whether to allow DNA-based descriptions (typifications) of species and by extension higher ranks of fungi, but what the precise requirements for such DNA-based typifications should be. We submit a tentative list of such criteria for further discussion. The present authors hope for a revitalized and deepened discussion on DNA-based typification, because to us it seems harmful and counter-productive to intentionally deny the overwhelming majority of extant fungi a formal standing under the International Code of Nomenclature for algae, fungi, and plants.

Effects of nitrogen deposition on carbon and nutrient cycling along a natural soil acidity gradient as revealed by metagenomics.

Nitrogen (N) deposition and soil acidification are environmental challenges affecting ecosystem functioning, health, and biodiversity, but their effects on functional genes are poorly understood. Here, we utilized metabarcoding and metagenomics to investigate the responses of soil functional genes to N deposition along a natural soil pH gradient. Soil N content was uncorrelated with pH, enabling us to investigate their effects separately. Soil acidity strongly and negatively affected the relative abundances of most cluster of orthologous gene categories of the metabolism supercategory. Similarly, soil acidity negatively affected the diversity of functional genes related to carbon and N but not phosphorus cycling. Multivariate analyses showed that soil pH was the most important factor affecting microbial and functional gene composition, while the effects of N deposition were less important. Relative abundance of KEGG functional modules related to different parts of the studied cycles showed variable responses to soil acidity and N deposition. Furthermore, our results suggested that the diversity-function relationship reported for other organisms also applies to soil microbiomes. Since N deposition and soil pH affected microbial taxonomic and functional composition to a different extent, we conclude that N deposition effects might be primarily mediated through soil acidification in forest ecosystems.

Global diversity and distribution of nitrogen-fixing bacteria in the soil.

Our knowledge of microbial biogeography has advanced in recent years, yet we lack knowledge of the global diversity of some important functional groups. Here, we used environmental DNA from 327 globally collected soil samples to investigate the biodiversity patterns of nitrogen-fixing bacteria by focusing on the nifH gene but also amplifying the general prokaryotic 16S SSU region. Globally, N-fixing prokaryotic communities are driven mainly by climatic conditions, with most groups being positively correlated with stable hot or seasonally humid climates. Among soil parameters, pH, but also soil N content were most often shown to correlate with the diversity of N-fixer groups. However, specific groups of N-fixing prokaryotes show contrasting responses to the same variables, notably in Cyanobacteria that were negatively correlated with stable hot climates, and showed a U-shaped correlation with soil pH, contrary to other N-fixers. Also, the non-N-fixing prokaryotic community composition was differentially correlated with the diversity and abundance of N-fixer groups, showing the often-neglected impact of biotic interactions among bacteria.

Global patterns in endemicity and vulnerability of soil fungi.

Fungi are highly diverse organisms, which provide multiple ecosystem services. However, compared with charismatic animals and plants, the distribution patterns and conservation needs of fungi have been little explored. Here, we examined endemicity patterns, global change vulnerability and conservation priority areas for functional groups of soil fungi based on six global surveys using a high-resolution, long-read metabarcoding approach. We found that the endemicity of all fungi and most functional groups peaks in tropical habitats, including Amazonia, Yucatan, West-Central Africa, Sri Lanka, and New Caledonia, with a negligible island effect compared with plants and animals. We also found that fungi are predominantly vulnerable to drought, heat and land-cover change, particularly in dry tropical regions with high human population density. Fungal conservation areas of highest priority include herbaceous wetlands, tropical forests, and woodlands. We stress that more attention should be focused on the conservation of fungi, especially root symbiotic arbuscular mycorrhizal and ectomycorrhizal fungi in tropical regions as well as unicellular early-diverging groups and macrofungi in general. Given the low overlap between the endemicity of fungi and macroorganisms, but high conservation needs in both groups, detailed analyses on distribution and conservation requirements are warranted for other microorganisms and soil organisms.

Long-term effects of forest fires on fungal community and soil properties along a hemiboreal Scots pine forest fire chronosequence.

We studied long-term effects of forest fires on the dynamics of soil fungal community along a post-fire chronosequence in hemiboreal Scots pine stands in north-western Estonia. Effects of fire on soil and fungi were studied on six sites that differed in time since fire (10, 21, 36, 67, 78 and 181 years ago), without further management interventions. Soil fungal communities along the chronosequence were dominated by soil saprotrophs and ectomycorrhizal (EcM) fungi. Across the chronosequence, the most dominant phylum was Ascomycota. The most abundant OTUs were identified as Umbelopsis sp., Hyaloscyphaceae sp. and Pezoloma ericae with relative abundances of 9.5, 8.9 and 6.8 %, respectively. Fungal species richness was similar among sample areas except in the area where fire occurred 36 years ago, where it was significantly lower. There were considerable differences in EcM fungal species composition along the chronosequence. The most recently burned site had Piloderma sphaerosporum, Pseudotomentella sp. and Clavulinaceae sp. as most abundant EcM OTUs while in three oldest burned areas Clavulinaceae sp. and Cortinarius sp. were abundant. Soil C and N stocks were lower in the most recently burned area but differences with other areas were not statistically significant. Soil pH had a significant effect on fungal species composition. Older areas had substantially lower pH compared to more recently burned areas.

Best practices in metabarcoding of fungi: From experimental design to results.

The development of high-throughput sequencing (HTS) technologies has greatly improved our capacity to identify fungi and unveil their ecological roles across a variety of ecosystems. Here we provide an overview of current best practices in metabarcoding analysis of fungal communities, from experimental design through molecular and computational analyses. By reanalysing published data sets, we demonstrate that operational taxonomic units (OTUs) outperform amplified sequence variants (ASVs) in recovering fungal diversity, a finding that is particularly evident for long markers. Additionally, analysis of the full-length ITS region allows more accurate taxonomic placement of fungi and other eukaryotes compared to the ITS2 subregion. Finally, we show that specific methods for compositional data analyses provide more reliable estimates of shifts in community structure. We conclude that metabarcoding analyses of fungi are especially promising for integrating fungi into the full microbiome and broader ecosystem functioning context, recovery of novel fungal lineages and ancient organisms as well as barcoding of old specimens including type material.

Structure and function of the soil microbiome underlying N2O emissions from global wetlands.

Wetland soils are the greatest source of nitrous oxide (N2O), a critical greenhouse gas and ozone depleter released by microbes. Yet, microbial players and processes underlying the N2O emissions from wetland soils are poorly understood. Using in situ N2O measurements and by determining the structure and potential functional of microbial communities in 645 wetland soil samples globally, we examined the potential role of archaea, bacteria, and fungi in nitrogen (N) cycling and N2O emissions. We show that N2O emissions are higher in drained and warm wetland soils, and are correlated with functional diversity of microbes. We further provide evidence that despite their much lower abundance compared to bacteria, nitrifying archaeal abundance is a key factor explaining N2O emissions from wetland soils globally. Our data suggest that ongoing global warming and intensifying environmental change may boost archaeal nitrifiers, collectively transforming wetland soils to a greater source of N2O.

Highly comparable metabarcoding results from MGI-Tech and Illumina sequencing platforms.

With the developments in DNA nanoball sequencing technologies and the emergence of new platforms, there is an increasing interest in their performance in comparison with the widely used sequencing-by-synthesis methods. Here, we test the consistency of metabarcoding results from DNBSEQ-G400RS (DNA nanoball sequencing platform by MGI-Tech) and NovaSeq 6000 (sequencing-by-synthesis platform by Illumina) platforms using technical replicates of DNA libraries that consist of COI gene amplicons from 120 soil DNA samples. By subjecting raw sequencing data from both platforms to a uniform bioinformatics processing, we found that the proportion of high-quality reads passing through the filtering steps was similar in both datasets. Per-sample operational taxonomic unit (OTU) and amplicon sequence variant (ASV) richness patterns were highly correlated, but sequencing data from DNBSEQ-G400RS harbored a higher number of OTUs. This may be related to the lower dominance of most common OTUs in DNBSEQ data set (thus revealing higher richness by detecting rare taxa) and/or to a lower effective read quality leading to generation of spurious OTUs. However, there was no statistical difference in the ASV and post-clustered ASV richness between platforms, suggesting that additional denoising step in the ASV workflow had effectively removed the 'noisy' reads. Both OTU-based and ASV-based composition were strongly correlated between the sequencing platforms, with essentially interchangeable results. Therefore, we conclude that DNBSEQ-G400RS and NovaSeq 6000 are both equally efficient high-throughput sequencing platforms to be utilized in studies aiming to apply the metabarcoding approach, but the main benefit of the former is related to lower sequencing cost.

Surveying for Batrachochytrium dendrobatidis and B. salamandrivorans in wild and captive amphibian populations in Estonia and Latvia.

The pathogenic chytrid fungi Batrachochytrium dendrobatidis (Bd) and B. salamandrivorans (Bsal) cause infections that have become primary drivers of amphibian biodiversity loss. While globally widespread, the distribution margins of Bd and Bsal have not been determined, and the presence of these pathogens has probably gone unnoticed in many areas, especially in northern Eurasia. To better understand the presence and distribution of both pathogens in the northern temperate and boreal forest biomes, 243 individuals were sampled from 8 native amphibian species across Estonia. Additionally, 68 amphibians were sampled from captive collections in Estonia and Latvia. Pathogen infection was assessed using metabarcoding of the ITS2 marker. No positive matches for Bsal infection were found. Bd was detected in 13 specimens, 3 of which were sampled at the Riga Zoo (with a prevalence of 5.2%) and 10 in natural environments in Estonia (3.3%). The infected wild individuals belonged to 6 amphibian species and were detected throughout the mainland of Estonia, but not on islands. Prevalence of infection with Bd ranged between 3.1 and 12.5% among native species. In addition, we found molecular evidence for a potentially new sister species to Bd in nature. Although outbreaks of chytridiomycosis have never been observed in Estonia, it cannot be excluded that the dynamics of local amphibian populations are affected by Bd infections. Therefore, further work, including capture-mark-recapture studies and long-term monitoring, are required to clarify the impact of Bd on amphibians in Northern Europe.

Perspectives and Benefits of High-Throughput Long-Read Sequencing in Microbial Ecology.

Short-read, high-throughput sequencing (HTS) methods have yielded numerous important insights into microbial ecology and function. Yet, in many instances short-read HTS techniques are suboptimal, for example, by providing insufficient phylogenetic resolution or low integrity of assembled genomes. Single-molecule and synthetic long-read (SLR) HTS methods have successfully ameliorated these limitations. In addition, nanopore sequencing has generated a number of unique analysis opportunities, such as rapid molecular diagnostics and direct RNA sequencing, and both Pacific Biosciences (PacBio) and nanopore sequencing support detection of epigenetic modifications. Although initially suffering from relatively low sequence quality, recent advances have greatly improved the accuracy of long-read sequencing technologies. In spite of great technological progress in recent years, the long-read HTS methods (PacBio and nanopore sequencing) are still relatively costly, require large amounts of high-quality starting material, and commonly need specific solutions in various analysis steps. Despite these challenges, long-read sequencing technologies offer high-quality, cutting-edge alternatives for testing hypotheses about microbiome structure and functioning as well as assembly of eukaryote genomes from complex environmental DNA samples.

Diversity and substrate-specificity of green algae and other micro-eukaryotes colonizing amphibian clutches in Germany, revealed by DNA metabarcoding.

Amphibian clutches are colonized by diverse but poorly studied communities of micro-organisms. One of the most noted ones is the unicellular green alga, Oophila amblystomatis, but the occurrence and role of other micro-organisms in the capsular chamber surrounding amphibian clutches have remained largely unstudied. Here, we undertook a multi-marker DNA metabarcoding study to characterize the community of algae and other micro-eukaryotes associated with agile frog (Rana dalmatina) clutches. Samplings were performed at three small ponds in Germany, from four substrates: water, sediment, tree leaves from the bottom of the pond, and R. dalmatina clutches. Sampling substrate strongly determined the community compositions of algae and other micro-eukaryotes. Therefore, as expected, the frog clutch-associated communities formed clearly distinct clusters. Clutch-associated communities in our study were structured by a plethora of not only green algae, but also diatoms and other ochrophytes. The most abundant operational taxonomic units (OTUs) in clutch samples were taxa from Chlamydomonas, Oophila, but also from Nitzschia and other ochrophytes. Sequences of Oophila "Clade B" were found exclusively in clutches. Based on additional phylogenetic analyses of 18S rDNA and of a matrix of 18 nuclear genes derived from transcriptomes, we confirmed in our samples the existence of two distinct clades of green algae assigned to Oophila in past studies. We hypothesize that "Clade B" algae correspond to the true Oophila, whereas "Clade A" algae are a series of Chlorococcum species that, along with other green algae, ochrophytes and protists, colonize amphibian clutches opportunistically and are often cultured from clutch samples due to their robust growth performance. The clutch-associated communities were subject to filtering by sampling location, suggesting that the taxa colonizing amphibian clutches can drastically differ depending on environmental conditions.

Tree Species Richness and Neighborhood Effects on Ectomycorrhizal Fungal Richness and Community Structure in Boreal Forest.

Tree species identity is one of the key factors driving ectomycorrhizal (EcM) fungal richness and community composition in boreal and temperate forest ecosystems, but little is known about the influence of tree species combinations and their neighborhood effects on EcM communities. To advance our understanding of host plant effects on EcM fungi, the roots of silver birch, Scots pine, and Norway spruce were analyzed using high-throughput sequencing across mature boreal forest exploratory plots of monocultures and two- and three-species mixtures in Finland. Our analyses revealed that tree species identity was an important determinant of EcM fungal community composition, but tree species richness had no significant influence on EcM fungal richness and community composition. We found that EcM fungal community composition associated with spruce depends on neighboring tree species. Our study suggests that at a regional-scale tree species identity is the primary factor determining community composition of root-associated EcM fungi alongside with tree species composition effects on EcM fungal community of spruce in mixed stands.

Temperature and pH define the realised niche space of arbuscular mycorrhizal fungi.

The arbuscular mycorrhizal (AM) fungi are a globally distributed group of soil organisms that play critical roles in ecosystem function. However, the ecological niches of individual AM fungal taxa are poorly understood. We collected > 300 soil samples from natural ecosystems worldwide and modelled the realised niches of AM fungal virtual taxa (VT; approximately species-level phylogroups). We found that environmental and spatial variables jointly explained VT distribution worldwide, with temperature and pH being the most important abiotic drivers, and spatial effects generally occurring at local to regional scales. While dispersal limitation could explain some variation in VT distribution, VT relative abundance was almost exclusively driven by environmental variables. Several environmental and spatial effects on VT distribution and relative abundance were correlated with phylogeny, indicating that closely related VT exhibit similar niche optima and widths. Major clades within the Glomeraceae exhibited distinct niche optima, Acaulosporaceae generally had niche optima in low pH and low temperature conditions, and Gigasporaceae generally had niche optima in high precipitation conditions. Identification of the realised niche space occupied by individual and phylogenetic groups of soil microbial taxa provides a basis for building detailed hypotheses about how soil communities respond to gradients and manipulation in ecosystems worldwide.